Strangelet spectra from type II supernovae 
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We study in this work the fate of strangelets injected as a contamination in the tail of a 
"strange matter-driven" supernova shock. A simple model for the fragmentation and braking of 
the strangelets when they pass through the expanding oxygen shell is presented and solved to un- 
derstand the reprocessing of this component. We find that the escaping spectrum is a scaled-down 
version of the one injected at the base of the oxygen shell. The supernova source is likely to produce 
low-energy particles of A ~ 100 — 1000 quite independently of the initial conditions. However, it 
is difficult that ultrarrelativistic strangelets (such as the hypothetical Centauro primaries) can have 
an origin in those explosive events. 
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I. INTRODUCTION 

More than a decade ago a celebrated paper by Wit- 
ten [Jj] (see also the previous works by Bodmer [g] and 
Chin and Kerman (3|) suggested that we may have been 
overlooking the true ground state of hadronic interac- 
tions. According to this work, strange matter (cold catal- 
ysed quark matter nearly symmetric in u 1 d and s flavors) 
would offer a possible form of bypassing the limitations 
imposed by the Pauli principle in ordinary matter be- 
cause of the existence of a third (strange) Fermi sea to 
share the energy of the system. Then, the energy per 
baryonic number unit would be smaller when compared 
to two-flavor quark matter and if this reduction is large 
enough, strange matter created after a weak-interaction 
time-scale may be the lowest energy state. Following 
this suggestion, which was essentially based on a bulk 
Fermi gas picture, the properties of strange matter and 
the droplet version, termed strangelets, were investigated 
Particular attention has been paid to the possible 
shell structure in the few-quark strangelets which 
are the ones expected to show up in heavy-ion collision 
experiments Q . A great deal of papers also explored as- 
trophysical scenarios which could render a non-zero ISM 
abundance of strangelets (which would cause all neutron 
stars to become strange stars) Jjjpl]-p^[ , the latter point 
being specially important because of the criticisms [ p~5|j3~6| 
raised to the SQM hypothesis based on pulsar glitch ob- 
servations, to which the basic strange star models [17] 
have no reasonable explanation to offer (see Refs. [IS] 
for more elaborated stellar models possibly containing 
an explanation of that data). 

Concerning the experimental detection at the Earth, 
it is interesting to note that, even before the official "in- 
vention" of SQM by Witten, some cosmic ray events |2C| ] 
were tentatively associated with quark blobs primaries 
pl| . A summary of the reported SQM candidates in sev- 
eral experiments is shown in Table 1. 

In Refs. QHIllHII as P ects of S Q M production 
have been investigated and discussed. One particularly 
puzzling aspect of all candidates is their relatively low 
baryonic number A. Calculations indicate [|] || that 
SQM tends to be more tightly bound for increasing A, 
so that injection of favored SQM fragments > 10 5 amu 
in any astrophysical event would require substantial re- 
processing to get down to ~ 10 2 amu strangelets. The 
authors of Ref . |26| have estimated the reprocessing time- 
scales by using model spallation cross sections of frag- 
ments of A amu with H and O of the form 
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where Ao is the preferred emitted cluster; and found that 
reprocessing is very ineffective, unless the strangelets can 
pass through a very dense oxygen shell. Although this sit- 
uation is very unlikely for the popular double-degenerate 
coalescence scenario jl 12|; it is precisely the situation 



expected in strange matter-driven supernovae scenarios 
p7| , p8| . In the latter a second shock arises because of the 
exothermic transition n — > uds + energy and is expected 
to carry a contamination of strangelets in the low- velocity 
tail as a byproduct of turbulent mixing |y],^8). Ejection 
of the strangelets with v ~ 0.1c is one of the possible 
mechanisms for a non-zero ISM abundance and thus com- 
petes with the coalescence events. The relative frequen- 
cies of both phenomena ~ 10 _2 yr _1 and ~ 10~ 4 yr -1 
respectively would be enough to identify the dominating 
source if the total mass ejected in strangelets could be 
calculated. However, since we lack of reliable estimates 
of the latter we have tried to infer the mass working back- 
wards from the reported events [14]. We have obtained 
1(T 6 M Q and 3 x KT 13 M by normalizing the flux to 
the events reported in Ref. |23| and Centauros (2(J re- 
spectively. However, these results refer strictly to the 
abundance of SQM primaries with A ~ 100 measured at 
the Earth, and therefore the question of the most likely 
" isotope" escaping from whatever source to the ISM still 
remains. 



II. SPALLATION OF STRANGELETS BY 
OXYGEN 

Let us address the specific case of SQM supernova ejec- 
tion. As explained above, the basic picture postulates 
of a (mildly relativistic to non-relativistic) strangelet 
gas travelling in the tail of the secondary shock, which 
encounters the expanding oxygen shell. According to 
Ref. p9[ , the expansion of the dense oxygen may be 
modelled by a number density evolution of the form 
n(t) — no exp(— £ / T exp ) where T exp ~ 0.1 s is a multiple 
of the free-fall time-scale appropriate for that physical 
situation. The (heavy) strangelets encountering oxygen 
targets will loose energy and baryon number, a process 
that may be described in the hydrodynamical approxi- 
mation by the following set of coupled equations 
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where v is the velocity of the strangelets relative to the 
expanding oxygen shell, m ox is the mass of an oxygen 
nucleus, Ro is the reference radius corresponding to an 
.4 = 200 strangelet (taken to be 5.8 fm) and C is the 
von Karman constant. We note that the r.h.s. of eq.(3) 
has been integrated over the fragment distribution and 
that we have included in the parameter < £ < 1 
the details of the transfer in the reaction A + 16 O — > 
A' + 16 0' + A + energy. 

Note that we have neglected the absorption of oxygen 
nuclei by the strangelet. These fusion reactions can be 
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crudely described by a geometric cross-section times the 
well-known Gamow factor. A simple calculation shows 
that the quotient of the fusion to the spallation cross- 
sections is ~ exp(— 327ra^4 1 / 3 cu~ 1 ), where a is the fine 
structure constant and the approximation Z ~ 2A 1 / 3 
has been made fl9|| . Absorption is thus suppressed by a 
large factor and is never important in the process we are 
considering. 

Eqs.(||) and (||) can be combined into a single differ- 
ential equation for the momentum loss which has the 
solution 
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where Vi and mi are the initial velocity and mass of the 
strangelet and D = (irC m ox i?g)/(2A (Jo) + £. Going 
back to eqs. (|^) and (^) we obtain the evolution of the 
strangelet mass with time as 
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with Ti = (nQViCr) 1 the initial mean-free path of 
the strangelets in the oxygen shell and fit) = (1 — 

exp(-t/Texp)). 

Finally, we find using eqs.(Q) and (g) the kinetic energy 
of the strangelet 
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and therefore the total energy of the strangelet Ej- from 
E$ = m 2 + E\, 



follows immediately that the strangelets can never sat- 
isfy the mass freezeout condition imposed on E^™' be- 
cause there is always enough energy available to power 
the spallations. Thus, the particles in this regime evapo- 
rate completely and do not constitute an astrophysically 
interesting case. 



B. 1/2 < D < 1 

In this case the kinetic energy decreases although the 
velocity increases with time. However, in the center-of- 
mass frame the kinetic energy E c ^ m - scales approximately 
as 
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(with E i ox the initial kinetic energy of the oxygen) 
and can not decrease either. Therefore the particles also 
evaporate in this regime. 



C. 1 < D < 4/3 

Now both the kinetic energy and the velocity of the 
strangelets decrease with time. Spallation proceeds until 
the point in which (see eq.7) the freezout condition is 
reached 
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This regime gives rise to a scaling law for the mass of 
the form 



III. REGIMES OF FRAGMENTATION 

An inspection of eqs. (4), (5) and (6) reveals that there 
are different regimes of fragmentation separated by the 
value of the exponent D. In all the four cases to be 
discussed below the mass of the strangelet decreases (i.e. 
spallation occurs) irrespective of the specific value of D; 
until the available kinetic energy in the center-of-mass 
frame E^ m ' is not enough to strip a chunk whick is bound 
by an amount Ao Eb/m p (Eb ~ WMeV is the binding 
energy per baryon number unit in the strangelet). The 
strangelet will survive as long as the latter condition can 
be reached before the mass is driven to zero (or, more 
precisely, to a value below m m j„ corresponding to the 
minimum stable baryon number of the strange matter 
A m in)- The possible regimes are 
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trip being the final or freezeout mass of the surviving 
strangelet. 



D. D > 4/3 

The situation is essentially the same as in the former 
point, but now the mass decreases slowly than before 
(see eq.5). This does not lead to any substantial change 



in the picture because the timescale Ti 



10 
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s IS so 



short that the strangelets would not be able to travel 
large distances before freezeout sets in or the oxygen shell 
expands substantially. 



A. < D < 1/2 

If D belongs to this range, eqs. (4) and (6) show that 
both the kinetic energy and velocity grow with time. It 
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IV. DISCUSSION 
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From the above expressions the full emerging spectrum 
in mass and energy could be computed for a given injec- 
tion spectrum at the source. Even without performing a 
detailed computation some general features are apparent. 
First of all eqs. (|)-9) show the dependence of the results 
with the parameter D, which are dramatically different 
depending on its actual value. Recalling the definition 
and using C = 0.5 as measured in the case of subsonic 
spheres in laboratory we have found a strict lower bound 
of D > 0.75 by setting £ = 0. Thus, it is quite likely 
that the physical situation corresponds to the third and 
fourth cases above. An important corollary is that the 
escaping spectrum of strangelets is a scaled-down ver- 
sion of the one injected by the secondary shock. A sec- 
ond major point is that the stripped fragments may or 
may not decay depending on the (unknown) actual value 
of the minimum stable strangelet A m i n . If A m i n ~ 10 
these fragments are likely to survive further spallation 
because of their non-relativistic character. On the other 
hand, they would quickly decay into ordinary hadrons 
if A m i n ~ 100. In any case the spallation interactions 
of the ejected strangelets with oxygen are so frequent 
in the expanding shell that the final distribution will be 
concentrated around the lowest stable " isotope" and the 
near the minimum (spallation cutoff) energy, although it 
is likely that further braking occurs either in the shell 
or in the ISM. For example, encounters of the escaping 
strangelets with giant molecular clouds are quite frequent 
and may help to strip a few baryon number units. Fi- 
nally, we note that is difficult in this scenario to obtain 
relativistic strangelets such as the ones required to fit the 
primaries of the Centauro events. The natural outcome 
from supernovae would be non-relativistic primaries in 
the range A ~ 100 — 1000; which are easily obtained 
from these events if the strangelets injected initially by 
the shock have 10 5 — 10 s amu for v% < 0.1 c (see eq. ||). 
A more realistic treatment of the initial strangelet spec- 
trum requires a model of fragmentation of the (initially 
homogeneous) strange matter fluid into strangelets in a 
turbulent environment [p8| ; even though we expect the 
mass scaling law eq.(9) also to hold in this case. This 
subject will be discussed in a future work. 

It is very important to measure the actual flux of can- 
didates to connect it with the details of the source and, 
through simulations like the one in Ref. [H , to measure 
the total mass in the galaxy (for example, a recent ex- 
periment seems to have measured a much lower flux than 
the one determined in p^] ). Further experimental and 
theoretical studies would be helpful to address the rela- 
tivistic ejection of strangelets and the form of the spectra 
for each astrophysical scenario. Work on these subjects 
is in progress and will be reported elsewhere. 
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TABLE I. SQM event candidates 



z 


A 


E/nuc [GeV] 


Ref. 




75 


~ 1000 


[20] 


20-40 




> 2.22 


[22] 


14 


400 


0.45 


[23] 


46 


> 1000 




[24] 



